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ABSTRACT 

We investigate the statistics of gravitational lenses in flat, low-density cos- 
mological models with different cosmic equations of state uj. We compute the 
lensing probabilities as a function of image separation 9 using a lens population 
described by the mass function of Jenkins et al. and modeled as singular isother- 
mal spheres on galactic scales and as Navarro, Frenk & White halos on cluster 
scales. It is found that COBE-normalized models with uj > —0.4 produce too few 
arcsecond-scale lenses in comparison with the JVAS/CLASS radio survey, a result 
that is consistent with other observational constraints on uj. The wide-separation 
(9 > 4") lensing rate is a particularly sensitive probe of both uj and the halo mass 
concentration. The absence of these systems in the current JVAS/CLASS data 
excludes highly concentrated halos in u < —0.7 models. The constraints can be 
improved by ongoing and future lensing surveys of > 10 5 sources. 

Subject headings: cosmology: theory - observations - cosmological parameters - 
gravitational lensing 

1. Introduction 

A number of observational data sets strongly suggest that the cosmological energy den- 
sity includes a component that is not associated with matter, either baryonic or dark. For 
example, measurements of the acoustic peaks in the cosmic microwave background (CMB) 
power spectrum point to a nearly spatially-flat cosmology (Pryke et al. 2001; Netterfield 
et al. 2001 and references therein). However, the matter density Q m inferred from cluster 
baryon fractions (e.g., David, Jones & Forman 1995) and galaxy redshift surveys (Percival 
et al. 2001 and references therein) falls far short of the critical value. Type la supernovae 
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studies offer independent evidence of a negative-pressure component of the energy density 
(Riess et al. 1998, 2001; Perlmutter et al. 1999). This energy component is described by the 
equation of state p = ujp with uj < 0. The value of uj is —1 for the spatially homogeneous 
cosmological constant A, and uj can be larger than —1 for other types of fields such as the 
quintessence (e.g., Frieman et al. 1995; Turner & White 1997; Caldwell et al. 1998). The 
latter clusters spatially on large scales, thereby modifying both the matter density fluctua- 
tion power spectrum and the CMB anisotropy (Ma et al. 1999). Large-scale structure and 
supernova observations currently favor — 1 < uj < —0.6 (Perlmutter, Turner & White 1999; 
Wang et al. 2000). The constraint can be improved with ongoing and future surveys based 
on, for example, the classical method of measuring the luminosity distance or the differential 
volume element as a function of redshift (e.g., Newman & Davis 2000). 

In this Letter we focus on strong gravitational lensing and examine the effect of the 
equation of state uj on the probabilities of producing multiply lensed systems. The number of 
expected lenses as a function of image separation provides a potential means of constraining 
uj because it is determined by factors such as the angular diameter distances to the lens and 
the source, the lensing cross sections, and the number density of the lenses, each of which 
depends on uj. This type of study is timely in view of the completion of the Jodrell-VLA 
Astrometric Survey (JVAS; e.g., Patnaik et al. 1992; King et al. 1999) and the Cosmic Lens 
All-Sky Survey (CLASS; e.g., Myers et al. 1999, 2001), which together provide the largest 
uniformly selected sample of radio lens systems and have yielded 18 new lenses thus far. The 
Sloan Digital Sky Survey (SDSS) will further increase the source population by one to two 
orders of magnitude. 

Since lensing probes the mass and not the light distribution, we model the lenses as 
a population of dark matter halos with an improved version of the Press-Schechter (1974) 
mass function. For galaxy-size lenses, we follow the tradition in lensing studies and use the 
singular isothermal spheres (SIS) as the mass profile (e.g., Turner et al. 1984; Narayan & 
White 1988; Kochanek 1996). The SIS ensures flat rotation curves and is consistent with 
constraints on the inner mass profiles of elliptical galaxies (e.g., Rusin & Ma 2001; Rix et al. 
1997; Kochanek 1995; Cohn et al. 2001). For cluster-size lenses, however, the mass profile is 
mostly determined by the dark matter. For this, we take advantage of the recent progress in 
high resolution A^-body simulations and model the lenses with the phenomenological profile 
of Navarro et al. (NFW, 1997). This approach allows us to calculate both small-separation 
(galaxy-based) and wide-separation (dark matter-based) lensing phenomena concurrently. It 
also enables us to relate the lensing probabilities directly to cosmological parameters through 
the mass power spectrum that governs the number density of lenses, and the lensing cross 
sections. Approximating lenses with a mixture of SIS and dark matter profiles is supported 
by simple baryon cooling models (Keeton 1998; Kochanek & White 2001) and has been used 
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recently to study the link between strong lensing and distant supernovae (Porciani & Madau 

2000) , halo density profiles (Keeton & Madau 2001; Takahashi & Chiba 2001; Wyithe et al. 

2001) , and Q m and Q\ (Li & Ostriker 2001). These studies, however, have assumed u = — 1. 
As we will show, relaxing the uo = — 1 assumption can have a significant effect on the lensing 
probabilities, in particular at wide angular separations. 



2. Lensing Probabilities and JVAS/CLASS Data 

In this Letter we consider cosmological models with a present-day density parameter of 
Q m = 0.35 in matter (with 0.05 in baryons and 0.3 in cold dark matter (CDM)), fl u = 0.65 in 
dark energy, and Hubble parameter h = 0.65. For the dark energy component, we analyze the 
standard cosmological constant (w = —1) and models with w < —1/3. This set of parameters 
is chosen to lie well inside the concordance region given by large scale structure observations 
(see Wang et al. 2000). Taking into account the current observational uncertainties in Vl m 
and h changes our constraints on u by up to ~ 10% (see § 3). 

The probability for lensing with an image separation greater than 9 is given by (e.g., 
Turner et al. 1984; Narayan & White 1988) 

/f z " dr r°° dn 
dz s V(z s ) dz t — dM—(M,z l )a lcns (M,z l )B(M,z h z s ), (1) 

J I J Mrnin 

where n(M, z{) is the physical number density of dark halos at the lens redshift z h <r lens (M, zi) 
is the lensing cross section of a halo of mass M at zi, B is the magnification bias (see below), 
and r is the proper cosmological distance to the lens, with dr/dz = cif " 1 (l + z)~ 1 { t Tl m (l + 
z) 3 + Q w (l + z)^ 1 ^]^ 1 / 2 for Q m + 0,^ — 1. The parameter M min is the mass needed for a 
halo at zi to create an angular separation 6 between the outermost images of a source at z s . 
By integrating over zi from to z s , we obtain the probability of lensing a source at z s . The 
last integral is over the redshift distribution V(z s ) of the sources (see below). 

For the halo number density, we use the recently calibrated expression dn/d In M = 
0.315 (1 + z) 3 (p/M)(dln5^ s /dkiM) exp(-| In ^ s + 0.61| 3 8 ) by Jenkins et al. (2001), where 
p is the mean mass density of the universe and 5 rms is the rms mass density fluctuation. This 
formula gives an improved match to the halo abundance found in numerical simulations 
compared with the classic Press & Schechter (1974) formula, which tends to overestimate 
the abundance of typical M* halos and underestimate massive halos. In the calculation of 
<5rms, we use the COBE normalized mass power spectrum of Ma et al. (1999) for the w > — 1 
models. The mass is taken to be the virial mass M = 47rA t) j r pi?^ r /3, where R V i r is the virial 
radius within which the average density is A vir p. For the standard CDM model with Q m = 1, 
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A vir is given by the familiar value A vir = 18n 2 ps 178; for flat ACDM models, A vir can be 
approximated by A vir (187r 2 + 82a; — 39:r 2 )/f2 m (,2) with re = f2 m (z) — 1 (Bryan & Norman 
1998). Jenkins et al. has specifically stated that their formula gives better fits to dn/dM 
with A vir = 178 regardless of the cosmological model, so we followed this instruction. 

The lensing cross section cri cns in eq. (1) depends on the mass profile of the lenses. Since 
our interest is to predict the lensing probability over a wide range of image separations, we 
consider both galaxy-size halos that are responsible for lens systems of a few arcseconds, 
and cluster-size halos for larger 6. In the former case, we approximate the galaxy mass 
distribution as an SIS with a density profile p(r) = a 2 / \2itGr 2 ) , where a v is the 1-d velocity 
dispersion. An SIS lens produces an image separation of 26 e, where 6e = 4n(o- v /c) 2 Di s /D s 
is the Einstein radius, and the cross section is cr lcns = ii(6 E Di) 2 = l67r 3 (a v /c) 4 (DiDi s /D s ) 2 
(e.g., Schneider et al. 1992), where D s ,Di, and As are the angular diameter distances to 
the source, to the lens, and between the lens and the source, respectively. Using a v = 
(irG 3 M 2 A vir p/Q) 1 / 6 , we can then relate a\ ens to the halo mass which is needed for eq. (1). 

For lensing by cluster-size halos, we approximate the cluster mass distribution by the 
density profile determined from halos in numerical simulations (Navarro et al. 1997): p{r) = 
p5u(r/r s ), where u(x) = l/[x(l + x) 2 }. This profile is shallower than the SIS for the inner 
parts of a virialized halo but is steeper at large radii. Other than the virial mass, this profile 
is described by a concentration parameter c = R V i r /r s , where R V i r is the virial radius of 
the halo discussed earlier, and r s is a scale radius. For c, we take into account both the 
mass and redshift dependence and use the relation of Bullock et al. (2001): c(M,z) = 
9(1 + z) _1 (M/M*)- a13 , where M* = 1.5 x 10 13 /i _1 M o . The density amplitude 5 is related 
to c by 8 = A 1 ,j r c 3 [ln(l + c) — c/(l + c)] _1 /3. The lensing cross section for NFW halos 
is determined by the parameter k = p5r s /S cr , where S cr = (c 2 /4:iiG)(D s /DiD[ s ) is the 
critical surface mass density. An NFW halo will have multiple images if the source is within 
the radial caustic of angular size l3 ra d from the halo center. We compute I3 ra d as a function 
of halo mass, lens redshift, and source redshift by solving df3/d6 = 0, where f3 and 6 are the 
positions of the source and the image, respectively. The cross section for lensing by an NFW 
halo is then <7i ens = n(l3 ra( iDi) 2 . We also need to determine the halo mass M min in eq. (1). 
Here we use the fact that the angular separation of the outermost images is insensitive to the 
value of (3 (Schneider et al. 1992) and simplify the calculations by using (3 = for a perfectly 
aligned source-lens configuration. We assume a "cooling mass" of M c ss 1.5 ■ lO 13 /i _1 M , 
above which the lenses are assigned the NFW profile and below which the lenses are SIS. 
Uncertainties introduced by M c on the lensing probability are discussed in § 3. 

To compare the predicted lensing probabilities with observational results, we use the 
combined data from JVAS and CLASS, which offer the largest uniformly selected sample of 
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gravitational lens systems. These surveys have discovered 18 lenses among a sample of ~ 
12000 flat-spectrum radio sources (Browne & Myers 2000; Browne et al. 2001; Myers et al. 
2001). A robust statistical analysis requires that careful cuts be made to the above sample, 
and this will be discussed in detail by Browne et al. However, because preliminary estimates 
indicate that the lensing rate will not differ much from the ~ 1/600 value derived here, we 
will assume a sample of 18 lenses and 12000 sources in the present analysis. 

Raw optical depths must be corrected to account for the magnification bias, which leads 
to an over-representation of lensed sources in any flux-limited sample (e.g., Turner et al. 1984; 
Maoz & Rix 1993). Magnification bias enhances the lensing probability of sources in a bin of 
total flux density (S) by the factor B(S) = _1 (S I ) J rf/x/i~ 1 P(/i)0(5'//i), where (fi(S) is the 
source luminosity function and P(fi) is the distribution of total magnifications (/z = £V 
where the magnification of the ith image is /ij) produced by the lens. The sources probed by 
CLASS are well-represented by a power-law luminosity function, (p(S) = dn/dS oc S~ v , with 
i] ~ 2.1 (Rusin & Tegmark 2001). The bias thus reduces to a simple form that is independent 
of flux density: B =< >. An SIS lens produces total magnifications described by the 
probability distribution -P(/i) = 8yU~ 3 ; so if r\ = 2.1, Bsis = 4.76. The cross sections of all SIS 
lenses are enhanced by this factor. The situation is more complicated for the NFW profile as 
its lensing properties depend on /t , which in turn depends on the halo mass and the angular 
distances to the halo and the source. We compute numerically the magnifications produced 
by NFW halos for 0.1 < n < 10 and use this to tabulate B nfw (k ). 

For the sources in the JVAS/CLASS survey, the redshift distribution is still poorly 
understood but the mean redshift is estimated to be < z s >= 1.27 (Marlow et al. 2000). 
Since the lower flux source distribution of the JVAS/CLASS survey is indistinguishable 
from the complete, brighter source distribution of the Caltech-Jodrell Bank VLBI sample 
(Henstock et al. 1997), we assume the latter quasar distribution for V(z s ) in eq. (1). 

3. Results and Discussion 

Fig. 1 compares the JVAS/CLASS data with the predicted lensing probabilities P{> 9) 
calculated from eq. (1) for the ACDM model and models with oo = —2/3,-1/2,-1/3. The 
probability at 9 > 4" decreases rapidly as oo increases towards 0. The most important 
systematics that affect wide separation lensing are due to the halo concentration parameter 
c(M,z). In Fig. 2a we quantify the dependence of P(> 4") on both u and the coefficient 
of c(M,z), c* = c(M*,z = 0). It shows that models with larger u can tolerate a higher 
halo concentration due to the lower lensing rates in these models. JVAS/CLASS thus far 
has detected no lens systems with 9 > 4" (Phillips et al. 2000). This excludes the highly 
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concentrated halos in uj < —0.7 models as indicated by the la, 90%, and 95% confidence 
contours in Fig. 2a. We also plot the la confidence levels that would be imposed if one wide 
separation lens system were discovered. We interpret the result as indicating that higher 
large separation lensing rates would lead to more refined constraints in the uj — c(M) plane. 

We can also compare the total number of lenses predicted by the models with the 18 
systems found among ~ 12000 sources in JVAS/CLASS. We find the expected number of 
lenses with 9 > 0.3" (the approximate angular resolution of JVAS/CLASS) to be 20.7, 18.2, 
15.5 and 9.6 for the ACDM, uj = —2/3, —1/2, and —1/3 models shown in Fig. 1, respectively. 
Because the SIS lensing cross section is several orders of magnitude higher than that of the 
NFW, the expected number of lenses is not sensitive to the concentration c(M). Instead, 
it depends more strongly on the cooling mass M c because a larger M c allows more halos to 
be modeled as SIS. In Fig. 2b we quantify the dependence of the total lensing rate on M c 
and uj by plotting the relative likelihood curves L(uo) oc p(uo) Nl (l — p(uo)) N ~ Nl for detecting 
Ni = 18 lenses from N = 12000 sources, where p(u) is the model-predicted lensing rate for a 
given uj. The range of M c is similar to that discussed in Kochanek & White (2001), with the 
upper and lower limits close to the cooling masses used by Porciani & Madau (2000) and Li 
& Ostriker (2001). The 2a confidence level in Fig. 2b suggests that the constraint on the 
equation of state is not sensitive to moderate alterations to the cooling mass. Models with 
uj > —0.4 are disfavored in all three cases. 

Note that our calculations implicitly assume that each halo below M c harbors sufficient 
baryons to give an SIS profile for the total mass. It is likely, however, that lower mass halos 
be devoid of baryons and therefore retain a shallower profile (e.g., NFW). Because an SIS 
lens has a much larger cross section than an NFW lens of the same mass, modeling all low- 
mass halos as SIS would overestimate the total lensing rate. Constraints on uj derived under 
our assumptions are therefore conservative - if low-mass halos actually contribute smaller 
cross sections than our calculations predict, more negative uj would be needed to reproduce 
the observed lensing rate. 

Both large separation lensing probabilities and the total predicted number of lenses are 
affected by the cosmological parameters h and Q m . Our investigations show that varying 
these parameters within the la error limits given in Netterfield et al. (2001) leads to almost 
parallel shifts in the contour lines and likelihood curves in Fig. 2 by 0.075 or less in uj, 
therefore not affecting the generality of our results. 

As the figures show, the predicted lensing probability decreases as w changes from — 1 
towards 0. The effect is particularly strong for lenses with 9 > 4", where the probability 
decreases by more than two orders of magnitude when uj is varied from —1 to —1/3. The 
dependence on uj enters in two ways in our calculations: one through the mass power spec- 
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trum upon which the halo density dn/dM depends; the other through kinematic factors 
such as the angular diameter distances which are functions of uj. For the power spectrum, 
different uj gives different linear growth rate for the density field. Specifically, for a fixed 
Q m and fl w , the growth becomes slower as uj increases towards zero because the energy den- 
sity in quintessence-type of fields dominates over that in matter at earlier times, resulting 
in earlier cessation of the gravitational collapse (Ma et al. 1999). The amplitude of the 
CO BE- normalized power spectrum is generally lower for larger u, resulting in a smaller halo 
density and hence a lower lensing probability. Besides the power spectrum, the dependence 
of the lensing probabilities on uj also enters through the effective lensing volume element 
(DiDi s /D s ) 2 (dr/dz), where the first factor comes from the lensing cross section <r lcns and the 
second factor is the distance to the lens in eq. (1). These effects are illustrated in Fig. 3. 

It is interesting to compare gravitational lensing studied in this Letter with classical cos- 
mological tests for constraining the equation of state. Ongoing efforts such as high redshift 
supernova searches and deep galaxy surveys offer promising ways to constrain uj by determin- 
ing the luminosity distance di and the cosmological volume element dV/dzdQ (e.g., Newman 
& Davis 2000; Turner 2001). In Fig. 3 we compare the luminosity distance and the volume 
element with the probability P(> 4") for large separation lensing. As we can see, P(> 4") 
depends more strongly on uj than the other quantities. Gravitational lensing statistics from 
ongoing and future surveys with 10 5 — 10 6 sources (e.g., the SDSS; see Cooray & Huterer 
1999) may offer an independent probe of uj, provided that systematics associated with input 
source redshift distributions and halo concentration parameters can be minimized. 

This paper would not have been possible without the hard work of everyone in the 
JVAS/CLASS team. We thank Mike Turner for useful comments. C.-P. M. acknowledges 
support of an Alfred P. Sloan Foundation Fellowship, a Cottrell Scholars Award from the 
Research Corporation, a Penn Research Foundation Award, and NSF grant AST 9973461. 
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Fig. 1. — Predicted lensing probability with image separation > 6 in cosmo logical models 
with different equation of state u. The lenses follow the mass function of Jenkins et al. (2001) 
and are modeled as SIS for M < 1.5 X 1O 13 /i _1 M and as NFW halos with the concentration 
parameter c(M,z) = 9(1 + z)~ l (M / 'MJ" - 13 , where M* = 1.5 x lO 13 /^ 1 ]^© (Bullock et al. 
2001). The histogram represents the survey results of JVAS/CLASS. 
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Fig. 2. — a) Wide separation lensing probability P(> 4") as a function of the equation 
of state us and the halo concentration c*, where c(M,z) = c*(l + z)~ 1 (M/M*) -0 ' 13 . Thick 
contours are various confidence upper limits for JVAS/CLASS, which has detected zero 
lenses with 9 > 4". (Upper-left corner is excluded.) Dotted contours represent the la 
confidence levels that would be imposed if one large separation lens were discovered. Dashed 
contours are for P(> 4") = 10~ 4 and 10~ 5 . b) Likelihood curves for detecting 18 lenses from 
12000 sources (found in JVAS/CLASS) for three different cooling masses. The dotted lines 
represent the la and 2a confidence levels. Models with uj > —0.4 are disfavored at 2a. 
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Fig. 3. — Ratio of four quantities in an w > —1 model relative to the ACDM model: 
the lensing probability P{> 4") for wide separation systems, the effective lensing volume 
element (DiDi s / D s ) 2 dr/dz (at Z\ = 0.5 and z s = 1.5), the volume element dV/dzdQ, and 
the luminosity distance d^ (both at z — 1). 



